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Abstract 



A ET-matrix formalism is used to relate the amplitudes for the three reactions 
pd — > ^Her/, vr" ^H— > ^Her/ and pd — > ^Htt^. Free parameters are fitted to the 
available experimental data and an extrapolation below the r/^He threshold is 
made to see the origin of the r/^He threshold enhancement. The existence of 
a virtual — and not a quasi-bound - state finds support in the data. The K- 
matrix permits a discussion of ry-vr mixing. A mixing parameter of 0.010(5), 
i.e. a mixing angle 9 = 0.6(3)°, is extracted from a best fit to the very recent 
pd — > ^Hevr'' reaction data. 
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I. INTRODUCTION 



In this paper we concentrate on the few-body interactions of 7] mesons in three nucleon 
systems. These complement our knowledge on the ry-nucleon interaction and properties of 
the T] meson itself. Considerable experimental, phenomenological and theoretical work has 
been devoted to understand the rj^He system. Early SATURNE experiments found a large 
cross section for the pd rj^'lie reaction close to the rj threshold, [1]. One interesting feature 
of this process is an enhancement of the meson formation amplitude in the few MeV energy 
region close to the threshold. A similar effect was also noted in the studies of the pn — > drj 
reactions made at CELSIUS, [2]. Both of these reactions indicate the possibihty for virtual 
or quasi-bound states to be formed in these systems. It is expected that in the deuteron 
the state of interest is a virtual state, whereas in helium it may also be virtual but there 
a bound state is not ruled out , [3], [4], [5]. While the slopes of the formation amplitudes 
indicate the existence of such states, more detailed properties may be found only with an 
extrapolation of these amphtudes below the ry^He threshold. This is made possible with the 
pion production experiments pd — > "^Htt"'" and pd — >-^He7r'^ undertaken by COSY, [6], [7]. 
These results obtained at backward angles supplement the older SATURNE cross sections 
measured some distance away from the threshold, [8]. 

In this work we present a multiple-channel X-matrix analysis of rj^He formation. In 
addition to the reactions hsted above, we also include the data from the tt"*" ^H— > ^He rj 
process studied at Brookhavcn, [9]. Altogether the data comprise 17 measurements. The 
number of important i^T- matrix parameters, obtained with some minor theoretical input, 
amounts to 4. Unfortunately the system is not strongly constrained, and so new data are 
welcome. 

A parallel study is devoted to the effect of rj-n mixing in these systems. It has been 
suggested in Ref. [6] that such a mixing is enhanced by the existence of a r/^He bound 
state. This enhancement has been found in a subsequent experiment, [7]. These two related 
questions are discussed in this paper, which is organised as follows: 

• Section II defines the mixed states of our system. 

• Section III discusses the i^-matrix method applied to the four-channel two body sub- 
systems of the meson-3N system. A full description, in the isospin symmetric limit, 
within a zero range approximation requires at least five real K-matnx elements. Cou- 
plings to the open three-body and four-body channels induce phases. Since the present 
data is not sufficient for a full determination of all the relevant i^-matrix elements, a 
model of S'(1535) dominance is used to remove some ambiguities. 

• In section IV the comparison with the experimental data indicates an attractive ?7^He 
scattering length and the existence of a rj^Re virtual state. Also the question of rj-ir 
mixing close to the r/^He threshold is discussed. 

• In section V some conclusions are made. 
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II. MIXED STATES 

Let the isospin invariant states of r] and 7r° mesons be described by | fj) and | tt). Due to 
some isospin mixing interaction Hm these states mix into the physical states | r]) and | tt). 
The relation of these two equivalent sets of states is 

\ 7^) = N[\ fj) - e \ Tt)] and \7r) = N[\Tr) + e \ f])], (1) 

where ^ is a mixing parameter and N — + 9"^ is the normalisation factor. These two 

sets of states form a complete orthonormal basis with the relations 

(r^ I tt) = = (r/ I 7f) and {f] \ tt) ^ NO ^ -{tt \ rj). (2) 

The mixing parameter follows from interactions at the quark level. The latter are due to 
differences in the light quark masses and to electromagnetic effects. For calculations of the 
mixing angle we refer to the review in Ref. [10]. This angle is usually expressed as 

^ {f]\H^\ Tt) 

where E"^, i?^ are the meson energies. 

The transitions in the few-body systems may be analysed in terms of the scattering 
matrices T that lead to the physical r] and tt in the final states 

T{r]) ^{pd\T\ rfHe) and T{t:) ^{pd\T\ n^He). (4) 

On the other hand, when discussing the formation processes and final state rj interactions, 
that are supposed to be isospin conserving, it is more convenient to use the isospin basis 

T{f)) ^{pd\T\ ffHe) and T{t^) ^{pd\T\ Tt^He). (5) 

The simple relationship between these amplitudes is 

T{v) = nn){n I V) + T(7r)(7r I 1^) ^ Tin) - eT{7r) (6) 

r(7r) = T{fj) (ry I tt) + T{7t) {n \ n) ^ T{7t) + eT{f]) . (7) 

The physical T matrix elements are used to describe the experimental data. However, to 
parameterise these data we are going to use the isospin basis amplitudes for the scattering 
matrix T and for the reaction matrix K. 

III. i^-MATRIX ANALYSIS OF THE rz-^HE PRODUCTION PROCESS 

In this section, the S-wave scattering matrix T for the idealised four channel two-body 
system | pd), pHer^), |^He7r'^), I^Htt"^) is introduced. These channels are denoted by the suf- 
fices p, ?7, TT, TT"*" respectively. The important couplings to open few-body channels are at first 
forgotten. The scattering matrix and the basic interactions are described and parameterised 
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(3) 



in terms of a real A'-matrix. Next, the Heitler equation for the T-matrix is solved. This 
becomes a simple matrix equation 



Tj^k = Kj^k + i^iKj^iqiTi,k, (8) 

where j, k, I are the channel indices and qi is a diagonal matrix of the CM momenta in each 
channel. Here the region of interest spans from about 10 MeV below the ?7^He threshold to 
some 10 MeV above it. For these low energies the experimental data exists. Also, in this 

region the K.j,k elements are believed to be constant. The model used here supposes the 
interactions to conserve isospin. In this way the fC- matrix is calculated in-between the | f]) , 
I 7f) and I 7f+) states. Since the isospin symmetry for the tt"*" is supposedly not violated one 
has I TT"*") =1 TT"*") and, in addition, some simple symmetries relate the matrix elements in the 
I^Hctt'^) and pH7r+) channels. In this way the parameters needed to fix the 4x4 X-matrix 
are reduced from 10 to 6. These arc Kp^p, Kp^,^, Kp^^,, -R'tt,??, -^r;,?;, -f^7r,7r- In practice, only five of 
these matrix elements are needed, since the coupling to the entrance pd channel is very small 
and so one finds Kp^p to be irrelevant. Other elements are related by the isospin symmetry: 

■^p,-jr+ V^-Kp^jr, K-jr+^jj ~ \/2-^7r,r;) -^7r+,7r+ '^^■K,-!n -^7r+,7r \/2-^7r,7r- 

The solution of the four dimensional Eq.(8) may be brought to a typical form 
where the scattering length in the rj^He system A^^^ is given by 



and the transition length Ap^^ is 



1 iSq^^Kj^.,^ 
For the pion production amplitudes one obtains 



A — I ^3g,r-^p,7r-^7r,?? , \ 



where 



-^p,7r[l '^qrjKri^ril ~\~ iqr^Kp^j^K^^,^ (13) 

1 iSq-i^K-jY^^ 

and the isospin relationship between the 7f+ and n is satisfied. Equations (9 - 13) contain ex- 
pressions which change rapidly in the small qjj region. One such term involves the scattering 
length Ari^ri- A quasi-bound state, if it exists, is given by the condition 

1 - iqrjAr^,,, = 0, (14) 

which is to be satisfied by a complex momentum qj^. For a large Aj^^rj this momentum is 
close to the threshold and the factor 1/(1 — ig^^r,,,?) induces rapid energy dependence of the 
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formation amplitudes in this region. This has been found in the rj production experiment 
[1], [3], but one can expect a similar behaviour in other channels. However, in the pion 
production amplitudes another factor arises which tends to suppress such an effect. It is 
due to the A^^^ of Eq.(13). There the existence of a quasi-bound state is expressed as the 
dominance of Kr^^^i- Therefore, if X^,^ were to differ from ^4^^^ in a significant way, then one 
would expect a sizable energy dependence in the pion formation amplitude. However, as 
will be seen later, the indications are that /T^,^ ~ A^,,,. This question is discussed in the 
following section. 



A. Estimate of large and small terms in the K-matrix 



The real physical situation involves the coupling of the two body channels to the contin- 
uum spectrum of three-body NNN, Ndn and four-body NNNn channels. The coupling to 
these systems may be strong and it is not easy to calculate. On phenomenological grounds 
it requires additional terms in the i^-matrix, so that 



1 - i-f^cc^c 



-K. 



c,k-i 



(15) 



where the summation (an integration) over the continuum few-body channels c is to be 
performed. This equation induces complex contributions to the real two channel matrix 
elements. In this subsection we estimate the magnitude of these contributions. The phases 
that arise are calculated in terms of a model or left to a best fit determination. The fine 
tuning of the parameters is done in the next section. The large matrix elements are those 
related to the low energy ?7^He channel. These enter essentially in the form of final state 
interaction factors in the r^^He system. One can visualise the relationship between the K- 
matrix parametrisation and a model description of the meson formation via an expression 



T 

■J- t 



(16) 



Here, ip denotes a wave function in the initial channel i and f/j^^ is an operator responsible for 
the meson formation. The term in square brackets is the final state wave function expressed 
in terms of the rjHe scattering matrix T^^^ = Ajj^j^/{1 — ig^A^^^). Up to terms linear in the 
final momentum g^, one obtains from Eqs.(9) and (16) 



A — TL 



1 + 



A, 



R 



(17) 



The Ui^r] and Ui^j^/R are results of the integrations in Eq.(16). The radius R refiects the 
range of final state interactions and is expected to be close to the ^He radius. Formula (17) 

contains effects from all the channels characterised by high intermediate momenta i.e. all K- 
matrix elements other than K^ j^ or A^ ,,, which have already been specified explicitly. Now, 
we show that the other matrix elements are small. Some will be needed only to the leading 
order, whereas others may be simply neglected. To do this we resort to the experimental 
data. 

The SATURNE cross section, [1], for the pd — >■ ^He rj reaction is given by T{p^ff) of 
Eq. (9). These data, given in Fig. 1, fix | Ap^^ \ rather precisely to the value of 0.013 fm. 
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On the other hand, provided | A^,,, |~ 5/m, there is a whole region of ReAj^^ri and ImAri^r] 
values that are equally likely. For cxamplc.the modulus of Eq.(9) is invariant with respect 
to the sign change ReA^ ^j —ReA^^^. These two possibilities describe different physics. 
Large positive values of ReA^^^^ correspond to a virtual state, an analogue of the NN spin 
state at low energies and the singularity of the scattering matrix given by Eq.(14) is located 
in the third quadrant of the complex plane. The other option, a negative length, signifies 
a quasi-bound state, which is analogous to the dcutcron. One sees that the formation cross 
section can not distinguish between these two possibilities. 

The Brookhaven data, for the 7r+ ^H— > ^He r] process, [9], permits one to extract two 
small X-matrix elements. These data are described by the amplitude 

1 K^^r, 

^-^ _ iq^Xr^^^){l - iq^K^^^) - qr^q-^Kl,^' ^^^^ 

giving I Kj^-ri |~ 0.07 fm, a value that depends only slightly on the choice of Kn^j^. The 
coupling of "^Her] to •^He(^H)7r is rather weak, and an inspection of Eq.(lO) tells us that 
its contribution to ImA^^^ is quite small. Therefore, the Brookhaven experiment implies 
that the rfYle state decays mainly into three- or four-body systems. Now, both r] formation 
experiments permit a simple description of the two body channels. Since it is well established 
that low energy rjN physics is dominated by the S'(1535) resonance we extend this dominance 
to few body systems. Thus the channel coupling is given by 

U,,, = F{q,)^,{^-^)^F{qj), (19) 

where the 7^ couple the resonance to the meson-nucleon channels, (. . .) denotes a suitable 
average of the resonance propagator over the binding and recoil energies and the F{qj) are 
the form factors for the meson-^He (^H) systems. The latter are expected to be about unity 
in the 1] channel and small in the tt channels due to the high momenta involved. The ratio 
■jTr/jri niay be extracted from 5'(1535) decay. Next, Eq.(19) — when combined with the 
(77,77) data — yield | K^^,^ \Ki 0.001 fm, a negligible value. There is another consequence 
of Eq.(19), the phase of ?7^^^ is given by the phase of the 5'(1535) propagator. For low 
energy 7/^He scattering, the relevant energies in the meson-nucleon system fall well below 
the resonance. Therefore, the dominant mode of decay is closed and f/^r is almost real. The 
uncertainty in the relative sign in the coupling constants yfjl and ^/^J may be removed by 
the 5'(1535) state wave function, where the SU(3) couphng mechanism of Ref. [11] gives a 
positive sign for X^.r;- 

We are now ready to study the pd ^'Hn^ reaction, [7]. Given a small value of K^^^^, 
alongwith a negligible i^Tr.Tr, then Eqs. (12, 13) and the experimental data from ref. [7] yield 
a crude estimate of I A„7r l~ 0.00021 fm. 



IV. FITTING THE iC-MATRIX TO EXPERIMENTAL DATA 

The data used consists of 8 measurements of the ^Her^ cross section in the threshold 
region [1], 4 measurements of the ^B.tt~^ cross section [7], 4 measurements of the ^HeTr" cross 
section [7], and one result for the tt"*" ^He rj reaction [9]. We limit the available data 
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to the low energy rj region to stay with the approximation of a constant 7^-matrix. The 
first two reactions yield absolute values of the Kp^^i and Kp .,^ matrix elements. However, the 
relative phase of these has to be left to an experimental determination. We set 

i^p,. = ^exp(z^), (20) 

where uj =\ Kp^jj/Kp^^r |= 5.54(50) is well determined from the rj and tt"'" formation exper- 
iments. The phase ip and the rj-ir mixing angle 6 are free parameters. To elucidate the 
interference pattern in the equations of the previous section, the tt formation amplitudes 
are now presented in a simplified form. Forgetting an irrelevant overall phase, up to terms 
linear in small X^,?? one has 

^j=Tip, 7f+) = KpAeMO^) + T^^^^l (21) 



T{p, TfO) = Kp,,[eMii^) + 7_ ^ ^]- (22) 

In these equations the two fC-matrix parameters are real. A positive sign for K^^^jj is preferred 
by the 5'(1535) dominance and also by the best fit to the experimental data. 

Using the Minuit minimization package, an overall best fit search to the data yields: 
Kp^r, = 0.0115(9) fm, Kp^^ = 0.00207(3) fm, X^,^ = 0.067(11) fm, K^^^ = , 
K^^r, = 4.24(29) + i0.69(81) fm, A^,^ = 4.24(29)' + i0.72(81) fm, 

i; = 4.14(27) and 9 = 0.010(0.005) = 0.6(3)°. The small difference between Im Kr,,n and 
Im Ari^ri is due to the explicit inclusion of Kj^^ri- As discussed in the last section, the large 
error in Im A^^^ arises since the r] formation cross section is not restrictive on the values 
of Ajj^jj. The real part of the ^He 77 scattering length is seen to be large and positive — 
signalling the existence of a virtual-state in this system. 



A. ri-TT mixing 

The idea behind the detection of ri-n mixing at COSY was to exploit the ratio of the 
charged and neutral pion cross sections Rmix, [6]- According to isospin invariance, Rmix 
should equal 2. However, the mixing induces corrections such that 



1 T(7r+) 


2 




2 




2 






2 


1 T{7r^) 1 


2 




TT 


)+T(7r)(7r 


\^)\ 


|2 N2 


1 + eT{f])/T{7T) |2 



(23) 



There is an additional reason for studying the ratio Rmix, rather than the separate cross 
sections, since in this way most of the systematic errors are removed. Because of that we 
limit our discussion to the data obtained in one laboratory. The amplitudes obtained in 
the previous section and 9 — 0.010(5) [6 = 0.6(3)"] reproduce the trends in the measured 
values to a fair degree as is shown in Table I and Fig. 2. In the latter, as a crude attempt 
to estimate systematic errors in the two measured cross sections, the effect of doubling the 
experimental error bars is shown as a dashed line. This is seen to have a minor effect. The 
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data and overall fit from the model display a maximum just above the r^-threshold. However, 
in the data a minimum below the threshold is also indicated, but it is not reproduced here. 
For comparison, in Fig. 3, the pd — > ^Hei:^ cross section calculated with a mixing angle of 
9 — 0.010 is plotted against the experimental results [7]. 

The last step in this analysis has been a more detailed extraction of the tj-tt mixing 
parameter 6. Notice that Eqs.(21) and (22), which arc essential for that procedure, have a 
very simple structure at the threshold, = 0. This point determines the unknown phase ip. 
We find two basic solutions that differ by the sign of 9, where the negative sign is ruled out 
on physical grounds. The actual value of 9 is then extracted by the energy dependence of 
the ratio Rmix a-nd the energy independence in | T(p, tt"*") | below the threshold indicated in 
Table I. However, the best fit parameters are not well determined and their errors are large. 
The need for more precise and more numerous data is evident. 

The value of the mixing parameter obtained here is smaller than the 9 = 1.5(4)° extracted 
from the nd — > rjNN reaction, [12]. Theoretical calculations yield values in the region of 
(0.75 — 0.85°) [13] - [15], although angles twice as large have also been suggested, [16]. 

B. Discussion of corrections 

So far we have assumed that all the isospin violation is due to the r]-7r mixing. However, 
there are other sources of such a violation. These can enter in two ways. 

• The two reactions pd ^Htt"'" and pd — > ^He7r° need not be exactly in the ratio 2 to 1 
as is assumed in Eq. (23), since they contain explicit isospin violation effects such as 
the difference between the ^He and ^H nuclear wavefunctions. Coulomb interactions 
and different meson momenta. This effect is analysed in terms of an additional free 
parameter A, by assuming | T'(7f+) P / | T(7t^) p to be 2A. The best fit value of 
A = 1.04(5) is found to improve our Rmix ratio in Fig. 3. In this way, the mixing 
parameter is reduced to ^ = 0.007(5). However, this procedure uses a very limited 
data base. An extension of the data could possibly lead to a different value of A. 

On the other hand qualitative theoretical arguments in Ref. [17] do seem to suggest a 
value of A that is greater than unity (pa 1.10). One conclusion from this type of overall 
renormalisation is that the error in 9 could well be larger. 

• An additional isospin violation effect in the pd — > ^Herj reaction could also arise from 
p-cu mixing. However, since our approach is based on phenomenological X-matrix 
parameters such an effect would not change directly the present determination of 9. 
Presumably this would contribute to an overall normalisation correction and so is 
taken into account by the above A correction. 

The conclusion is that, cvcnthough these two effects could be at a 10% level, they do not 
lead to any dramatic effect at the rj- threshold and so are incorporated in the multiplicative 
factor A. 

In addition to the above corrections, it should be remembered that Eq. ( 23) is written 
down for S-waves under the assumption that these dominate in the backward scattering, 
whereas a more correct expression would involve the effect of higher partial waves. That 
possibility was incorporated by simply adding an additional "background" contribution as 
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a complex constant (c) to T(p, vr). However, the best fit procedure indicated that c was very 
small and so ruled out any significant contribution of this kind. 



V. CONCLUSION 

The i^— matrix formalism developed here is able to account for the structure seen at the 
77-threshold in the experimental ratio | T'(7f+) p / | T(7f°) p as a manifestation of rj-n mixing. 
It also enables an estimate of 0.010(5) to be made of the mixing angle. Unfortunately, at 
present this estimate has a large uncertainty, which could be significantly reduced by the 
removal of several uncertain systematic effects in the available experimental data. This 
clearly exposes the need for more precise data over the energy range covering the T^-threshold. 
Such data should be detailed within the range piab— 1-55 - 1.59 GeV. In addition, some data 
points further from the threshold would be very valuable in order to study the non-threshold 
value of I T{'K~^) ^ / | r(7f°) ^ and so tie down more precisely model parameters such as A. 
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TABLES 



TABLE I. The experimental and calculated amplitudes | T{tt'^) p [10~^/m^] for the tt"^ pro- 
duction. Other cohimns give ratios R =| T{7:~^) p / | T{7t^) p calculated with a mixing angle 
of ^? = 0.010. The experimental results are from COSY [7]. The first column gives the proton 
laboratory momentum in GeV. 



Plab 


1 T{7T+) IL, 


1 n^^) \lalc 


Rexp 


-^calc 


1.560 


47.4(3.9) 


49.7 


2.05(0.17) 


2.04 


1.570 


45.8(6.6) 


48.8 


1.84(0.27) 


2.08 


1.571 


47.5(2.2) 


46.0 


2.24(0.11) 


2.17 


1.590 


62.6(8.6) 


47.6 


2.57(0.27) 


2.11 



Figure Captions 

Fig.l The pd — > rj^Re data (the formation amplitude | / | = | T(p, rj) \ squared ) of Ref. [1] 
and the corresponding X-matrix fit. 

Fig. 2 The ratio R^ix =\ T{n+) ^ / | T(7r°) p calculated with the mixing angle 6 = 0.010. 
The experimental results are from COSY [7]. The dashed curve is the AT— matrix fit, where 
the error bars of the pd — > ^Htt"*" and pd — > ^He7r° data have been doubled in an attempt to 
simulate systematic errors in that data. These results include the correction factor A = 1.04 
discussed in Subsect. IV B. 

Fig. 3 The pd — > ^Hei:^ cross section calculated with the mixing angle = 0.010 and 
A = 1.04. The experimental results are from COSY [7]. 
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FIGURES 

FIG. 1. The pd rj^He data (the formation amphtude | / |=| T{p,rf) \ squared ) of Ref. 
and the corresponding if-matrix fit 
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FIG. 2. The ratio Rmix =\ T{tt+) ^ / | T{Tr^) p calculated with the mixing angle 9 = 0.010. 
The experimental results are from COSY [7]. The dashed curve is the ii'— matrix fit where the 
error bars of the pd — > ^Htt"'" and pd — > ^HeTT*^ data have been doubled in an attempt to simulate 
systematic errors in that data. These results include the correction factor A = 1.04 discussed in 
Subsect. IV B 
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FIG. 3. The pd ^Hen^ cross section calculated with the mixing angle 9 = 0.010 and A = 1.04. 
The experimental results are from COSY [7]. 
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